m 
o 
o 

(N 
< 



Nonlocal mixing of supercurrents in Josephson ballistic point contact 

Yu.A. Kolesnichenko, A.N. Omelyanchouk, S. N. Shevchenko 
B.I. Verkin Institute for Low Temperature Physics and Engineering, 

National 

Academy of Sciences of Ukraine, J^l Lenin Ave., 61103, Kharkov, 

Ukraine 
(February 2, 2008) 



o 
o 

Oh- 



o 
o 



> 

o 
m 
o 



I 

o 
o 



X 



We study coherent current states in the mesoscopic su- 
perconducting weak link simultaneously subjected to the or- 
der parameter phase difference on the contact and to the 
tangential to the junction interface superfluid velocity Vs in 
the banks. The Josephson current-phase relation Ij{4>) con- 
trolled by the external transport current 7t(vs) is obtained. 
At (j> close to TT the nonlocal nature of the Josephson phase- 
dependent current results in the appearance of two vortexlike 
states in the vicinity of the contact. 

The superfluid flow of Cooper pairs in superconductor 
is related to the space dependence of the phase x of the 
order parameter. In (quasi)homogeneous current state 
the supercurrent density j locally depends on the super- 
fluid velocity Vs = 2^Vx(r). Such a state is realized 
in narrow films or wires [1]. In the case when the phase 
X strongly varies in the scale of superconducting coher- 
ence length ^0 the relation between the current density 
j(r) and x(r) becomes nonlocal. This situation (opposite 
to the homogeneous current state) is realized in Joseph- 
son weak links (for review see [2]), e.g. in superconduct- 
ing point contacts - microconstrictions between two bulk 
superconductors (banks). The Josephson current nonlo- 
cally depends on x(r) and is determined (parameterized) 
by the total phase difference (j) across the weak link. The 
current-phase relation j(0) for ballistic point contact was 
obtained in [3]. The nonlocal nature of Josephson cur- 
rent in mesoscopic junctions was demonstrated by Heida 
et. al. [4] and studied in theoretical papers [5]. 

The Josephson weak link could be considered as a 
"mixer" of two superconducting macroscopic quantum 
states in the banks. The result of the mixing is the 
phase dependent current carrying state with current flow- 
ing from one bank to another. The properties of this 
state depend on the properties of the states of the banks. 
For example, in Josephson junction between unconven- 
tional (d-wave) superconductors the surface current, tan- 
gential to the contact interface, appears simultaneously 
with Josephson current (see e.g. [6]). 

In this paper we study coherent current states in the 
Josephson weak link between conventional superconduc- 
tors, whose banks are in the homogeneous current states. 
The questions raised in our consideration are the fol- 
lowing ones: How two superconducting current carrying 
states in the banks are coherently mixed by a mesoscopic 
Josephson junction, or in other words, what is a result of 
the interplay between transport current jri'^s) flowing 



parallel to the junction interface and nonlocal Josephson 
current j j (</>)? How the Josephson properties of the sys- 
tem are influenced by the external controlling transport 
current? We have found that the distribution of the cur- 
rent in a region of nonlocal mixing strongly depends on 
the global phase difference between banks and for cj) — tt 
contains the vortexlike states. The current-phase relation 
jj{(f>) at near tt essentially depends on the superfluid 
velocity in the banks Vg, in particular, the absolute value 
of the derivative djj/dcj) at = tt is suppressed by the 
transport current. 

We consider the Josephson weak link with direct con- 
ductivity - a microbridge between thin superconducting 
films. The bridge sizes, length L and width 2a, are as- 
sumed to be smaller than the coherence length ^o- In this 
case even for temperature T near the critical temperature 
Tc the local description based on the Ginzburg-Landau 
approach is not applicable. To describe the coherent cur- 
rent states in the system we use the quasiclassical Eilen- 
berger equations [7], which are valid for temperatures 
< T < Tc and for arbitrary relation between the con- 
tact size and coherence length fo- On the other hand, 
we assume that a and L are much larger than the Fermi 
wavelength Xp. The electron mean free path is supposed 
to be much larger than ^o- 

Suppose the homogeneous transport current It with 
a superfiuid velocity flows in the banks of the con- 
tact. The situation with controlled phase difference 
and preset current It may be realized if the microbridge 
is incorporated in a cylindrical thin film (Fig 1). Let 
the radius of the cylinder be less than London penetra- 
tion depth and larger than the coherence length. In this 
case the phase difference </> is governed by the external 
magnetic flux ^, (j) — and the external transport 

current It flowing along the cylinder is homogeneously 
distributed far from the microconstriction. 

The Eilenberger equations for the ^-integrated Green's 
functions have the form [7] : 

d 

vj^— G^(vF,r) + [uj% + A(vF,r),G,^(vF,r)] = 0, 

(1) 



dr 
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T3 is the Pauli matrix, A is the superconducting order pa- 
rameter, and Gi^{'VF,r) is the matrix Green's function. 



FIG. 1. Scheme of the reahzation of the microbridge with 
transport current in the banks and controlled phase difference. 

which depends on the electron velocity on the Fermi sur- 
face vp, the coordinate r, and the Matsubara frequency 
LO ~ {2n + l)7rT, with n being an integer number. 

The order parameter A is determined by the self- 
consistency equation 

A(r)=^Ar^(/„(v^,r))^^. (3) 

Solution of the matrix equation (1) together with Eq.(3) 
determines the current density j(r) in the system 

j(r) = -27rieiV(0)r^(vF5„(v^,r))^^ . (4) 

Here A is the BCS coupling constant, A^(0) is the density 
of states at the Fermi surface, is the averaging over 

directions of the velocity vp. 

If the film thickness w is much smaller than ; the spa- 
tial distributions of A(r) and j(r) depend only on coor- 
dinates in the plane of the film and the Eilenberger equa- 
tions (1) reduce to the two-dimensional ones. We solve 
these equations in the model of the microbridge as a slit 
in thin impenetrable partition (L = 0) at a; = between 
two half- planes x ^ (Fig. 2). The equations (1) for 
Green's function Gcj{vp,x,y) have to be supplemented 
by the continuity condition at the slit {x = 0, |y| < a) and 
by the condition of the specular reflection at the line 
{x = 0, \y\ > a). Far from the constriction the Green's 
functions must satisfy to the conditions, which describe 
the homogeneous current parallel to the axis. 

As it was shown in [3] in the zero approximation on the 
small parameter o/^o <C 1 for a self-consistent calculation 
of the superconducting current it is not necessary to find 
the spatial dependence A{x,y). In the same approxima- 
tion the superfluid velocity does not depend on coor- 
dinates. The spatial variation of A and Vg is essential 
at the distances p from the contact p < a. The Green's 



FIG. 2. Model of the microbridge as the slit in the thin 
insulating partition 

functions are varied at the distances of order and in 
the main approximation on the parameter a/^o ^ 1 they 
are defined by the values of A (v_, ) and v, in the banks of 
the contact. For A and v,, being c;onstants at each half- 
plane an analytical solution of Eilenberger equations can 
be found by the method of integration along quasiclas- 
sical trajectories. Under the condition a/^o ^ 1 such 
solution is self consistent. In any point p = {x,y) all bal- 
listic trajectories can be classified as transit trajectories 
(marked by "1" at Fig. 2), for which v^^- e oi{p) {a{p) 
being the angle at which the slit is seen from the point 
p) and non-transit trajectories, 'Vp ^ a{p), (marked by 
"2" at Fig. 2). For transit trajectories the Green's func- 
tions satisfy the boundary conditions in the both banks. 
For the non-transit trajectories they satisfy to specular 
reflection condition at the partition and the conditions 
in the left or right bank. Making use of the solution of 
Eilenberger equations, we obtain the following expression 
for the current density (4) at the slit: 

j(a; = 0, |y| < o, (/),v,,) = (5) 

ifi sin ^ — 770; cos ^ \ 
47r|e|jV(0KT^(vJm \ ' . l ) , 
ooX r?f2cosf -iwsmf /„ 

where f2 = v^wM-A^, ui = iv + ippVg, v = vp/vp is the 
unit vector, r] = sign{vx)- We should require ReQ. > 0, 
which fixes the sign of the square root to be sign(pi?v,,). 
Under the condition a/^o 1 the current density at the 
slit does not depend on the ?/— coordinate, and the total 
current through the contact (Josephson current) is equal 
to I J = 2awjx{x 0, \y\ < a, (^,Vs). 

For Vs = the component of the current (5) tangential 
to the contact jy = and for the Josephson current den- 
sity jj = jx we have the result obtained in the paper [3] . 
In general case Vg ^ the current (5) has both jj and jy 
components. The tangential current jy depends on the 
phase (j) and is not equal to the transport current density 
]t in the banks. In particular, at (j) near tt it goes in 
opposite direction to the external transport current (see 
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FIG. 3. Josephson current Ij versus phase 0. 
/o = law ■ Att |e| N{0)vfTc. 

below) . 

To describe the influence of the transport cur- 
rent in the banks on the Josephson current we in- 
troduce the dimensionless parameter q = VsPf/^o 
(Aq = A(T = Q,Vs = 0)) . The value of q is varied in the 
range < g < gc- The critical value qc corresponds 
to the critical current density in the homogeneous cur- 
rent state. At zero temperature qc = 1, and the gap A 
does not depend on q [8]. In Fig. 3 we plot the Joseph- 
son current at temperature T = O.lTc for different 
values of q. The presence of the tangential transport 
current in the banks suppresses the value of the critical 
Josephson current and essentially changes the derivative 
dlj/d(j) at (j> = IT. We emphasize, that the dependence 
of the Josephson current Ij (0) on q, which is shown in 
Fig. 3, does not relate to the suppression of the gap by 
the transport current, which is negligible for such low 
temperatures. 

The derivative dij/dcj) at ^ = tt determines the ki- 
netic inductance of Josephson jimction [2], which is rele- 
vant, e.g., for squid's operation. The expression for the 
derivative of the Josephson current Ij at <p = tt has the 
form 

# mvs V ttT / 

where function f{x) is plotted in Fig. 4. The derivative 
dlj/d<p{(j) = tt) is inversely proportional to tig at T <C 
PfVs and it is inversely proportional to T at T prVg- 
By using the Green's functions along transit and non- 
transit trajectories, calculated in the main approxima- 
tion on the small parameter a/^o, we can find the spa- 
tial distributions of the order parameter and the current 
density in the contact (see Ref. [3]). The numerically cal- 
culated current density distributions for different values 
of phase </> and the temperature T = O.lTc are shown in 
Figures 5 and 6. For small values of the phase difference 
(/) between banks, the current density j(p) is just the vec- 
tor sum of the homogeneous transport current density in 
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FIG. 4. The function f{x) from expression (6) 
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FIG. 5. Vector plot of the current density for (j) = n/2 and 
q = 0.5. Numbers maxk x and y axes in the units of the 
contact size a. 

the banks jriys) and the conventional Josephson current 
p, = 0) (Fig. 5). For near tt the constructive in- 
terference of supercurrents takes place. At (f) = tt there 
are no Josephson current, jj = 0, and the current is 
distributed in the way that there are two antisymmetric 
"vortices" close to the contact region (Fig. 6). Far from 
the constriction (at the distances p ^ o) the inter- 
ference current is spread out and the current density is 
equal to its value in the banks. 

Simple and transparent expression for the current den- 
sity distribution j(/9) can be found for temperatures close 
to the critical temperature (Tc — T <C Tc). At the dis- 
tances from the contact, which are less than the coher- 
ence length ^0 and arbitrary in comparison with the size 
a, j(p) takes the form: 

j(P,0,Vs) =Jj{p,(l>) +jT(vs) +jjT(p,0,Vs); (7) 
ij{P,(l>) = 2jc sin (vsi.gn(?;:r))^eo;(p) 5 

ijT{p,<p,Vs) = jck{l - cost/)) (vt;j,)^g„(^) , 



\\\\ 



f \ 



FIG. 6. Vector plot of the current density for = tt and 
q = 0.5. 



where 



3c iT,v,) = 



n\e\N{0)vF {T,Vs) 



(8) 



is a critical current density of the contact at T w Tc, 
k = {U<;{3)/tt^){vsPf/Tc). We detach explicitly the 
Josephson current jj{p,(f>), and the spatially homoge- 
neous (transport) current density jT(vs) produced by the 
superfluid velocity v^, and write the total current (7) as 
the sum of three components: jj, jr, and the rest - the 
"interference" current jjT- The macroscopic quantum 
interference takes place in the vicinity of the contact re- 
gion where both coherent current densities and 
jT(vs) exist. We emphasize, that at near tt at the slit 
the "interference" current jj^ is antiparallcl to j^. If the 
phase difference (f> = n, the current jjx = — 2jT. When 
there is no phase difference (at </> = 0), we have jjr = 0. 

In conclusion, we have investigated the coherent cur- 
rent states in the Josephson ballistic point contact simul- 
taneously subjected to the order parameter phase differ- 
ence 4> and to the tangential to the junction interface the 
superfluid velocity Vg in the banks. The current-phase 
relation Ij {(p) is shown to be controlled by the transport 
superconducting current /t(vs). Thus, varying /^(vs), 
the characteristics of the weak link, such as the shape of 
the Josephson current-phase relation and the value of the 
critical current, can be changed. The similar effect can 
be produced by the increasing of the temperature T of 
the system. But, as compared to the controlling by the 
transport supercurrcnt, the increasing of the temperature 
leads to additional thermal noise. 

Moreover, the current distribution pattern in the vicin- 
ity of the contact was obtained. The current pattern 
drastically depends on the external phase difference 0; in 
particular, at = tt the existence of two antisymmetric 
vortexlike current structures is predicted. Considering 
the current pattern, we have also demonstrated that the 
superposition of the supercurrents in the vicinity of the 
weak link is not just their vector sum. These results can 
be relevant in a wide range of problems, in which the 



current (and corresponding magnetic field) distribution 
in the vicinity of the weak link is important. 
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